
CENTER FOR SPACE RESEARCH 
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 






This work was supported in part by the National Aeronautics 
and Space Administration under Grant NGL 22-009-372. 


STELLAR WINDS DRIVEN 
BY ALFVEN WAVES 

by 

John W. Belcher and Stanislaw Olbert 
CSR TR-73-5 October, 1973 

Department of Physics and Center for Space Research 
Massachusetts Institute of Technology 
Cambridge, Massachusetts 02139 


c 





ABSTRACT 


We consider models of stellar winds in which the dynamic 
expansion of a corona is driven by Alfven waves propagating 
outward along radial magnetic field lines. In the presence of 
Alfven waves, a coronal expansion can exist for a broad range 
of reference conditions which would, in the absence of waves, 
lead to static configurations. Wind models in which the accel- 
eration mechanism is due to Alfven waves alone exhibit lower 
mass fluxes and higher energies per particle as compared to 
wind models in which the acceleration is due to thermal processes. 
For example, winds driven by Alfven waves exhibit streaming ve- 
locities at infinity which may vary between the escape velocity 
at the coronal base and the geometrical mean of the escape 
velocity and the speed of light. We derive upper and lower 
limits for the allowed energy fluxes and mass fluxes associated 
with these winds . 
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I . INTRODUCTION 

The properties of stellar winds which are thermally 
driven have been the subject of extensive theoretical investi- 
gation over the last fifteen years. Dynamical models of such 
winds have become increasingly sophisticated, including effects 
due to the two fluid nature of the plasma, magnetic fields and 
stellar rotation, the inhibition of thermal conductivity, the 
propagation and damping of hydromagnetic waves in the expanding 
solar corona, and many others. For comprehensive reviews of 
these topics, see the articles by Parker (1971) and Barnes (1973), 

or the book by Hundhausen (1972) . 

Recently, Belcher (1971, 1972), and, independently, Alaz- 

raki and Couturier (1971) have considered modifications of 
polytrope wind models due to the presence of undamped Alfven 
waves propagating outward along radial magnetic field lines. 

The interaction of the waves with the streaiming plasma produces 
an outward pressure gradient, analogous to that of a radiation 
pressure, which results in a radial acceleration of the wind. 

In this manner, undamped wave energy fluxes propagating out- 
ward into an expanding corona are completely transformed into 
enhanced streaming energy fluxes of the wind at large distances 
from the star. Parker (1965) was the first to suggest that 
undamped Alfven waves could affect the dynamics of the solar 
wind. However, detailed considerations of the problem were 
not undertaken until observational evidence suggested that 
Alfven waves generated close to the sun are in fact present 
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at 1 a.u. (Coleman 1967, Unti and Neugebauer 1968, Belcher 
et al . 1969, Belcher and Davis 1971), Since the original one 
fluid treatments , various authors have also considered the 
effects of Alfvenic wave pressures in two fluid models (Hollweg 
1973a) , the modifications due to non-WKB terms in one fluid 
models (Hollweg 1973b) , and the effects of finite and large 
amplitude Alfven waves in one fluid models (Whang 1973, Barnes 
and Hollweg 1973) . Hollweg (1972) has also considered possible 
generation mechanisms for these waves in the solar chromosphere. 

For obvious historical reasons, the initial treatments 
of wave pressures are primarily concerned with situations in 
which the Alfven wave energy flux across the base of the corona 
is less than the conductive flux of thermal energy - that is, 
the addition of wave pressures is considered to be a modification 
of an essentially thermal process. In the opposite extreme, 
however, Alazraki and Couturier (1971) point out that wind 
solutions to the equations of motion always exist as long as 
the wave energy flux is non-zero, even if the conductive flux 
of thermal energy is identically zero. Such wave driven winds 
may exhibit large energies per particle at infinity, in con- 
junction with small mass fluxes (Belcher 1971). It is thus 
possible for Alfven waves alone to drive a coronal expansion, 

^he properties of winds produced in this manner may be very 
different from those of thermally driven winds. In the present 
paper, we investigate in detail the characteristics of stellar 
winds which are primarily driven by low-frequency, outwardly- 
propagating Alfven waves generated close to a star. The thermal 
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properties of the plasma are represented by a polytrope relation 
between density and pressure, and we consider only thermal 
parameters which in the absence of waves would result in static 
atmospheres. For a given set of initial parameters {wave 
strengths, temperatures, densities, and so on) at some reference 
level close to the star, we wish to determine whether or not 
the dynamical expansion of the atmosphere into a stellar wind is 
possible, and, if so, to ascertain the mass and energy fluxes 
associated with that wind. We assume that all generation 
mechanisms for the Alfven waves (e.g. , convective zones) occur 
inside the reference level, and that there is no damping of the 
waves external to the reference level. As we shall see, there 
are a broad range of wind solutions possible. 

Before proceeding with the detailed mathematics, we offer 
some rationale for the formulation we use. First, to keep the 
calculation tractable, we consider only radial streaming in the 
presence of a radial magnetic field, with no stellar rotation. 
Second, the winds we shall encounter may exhibit rapid decreases 
in mass density outward from the reference level, whereas the 
magnetic field strength decreases less rapidly, as 1/r^. Con- 
sequently, the Alfven velocity in the low density, field 
dominated plasma may be high, and we must insure that it does 
not exceed c, the speed of light. As we shall show, the fact 
that both Alazraki and Couturier (1971) and Belcher (1971) 
allow the Alfven velocity to be arbitrarily large invalidates 
their results at low mass fluxes. Third, in some limits we 
shall encounter winds with arbitrarily large energies per particle 
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far from the star, in conjunction with very low mass fluxes. 
Although the physical validity of these solutions is question- 
able, we must allow for the possibility of relativistic streaming 
velocities at infinity to handle the limits properly. Finally 
we shall find circumstances in which wind solutions formally 
exist even for tightly bound atmospheres near massive objects 
with high escape velocities. Thus our initial approach should 
allow for escape velocities, Alfven velocities, and radial 
streaming velocities which may be comparable to the speed of 
light. We consider only situations in which the local sound 
velocity and the transverse velocity perturbation associated 
with the wave are small compared to the speed of light. To 
insure the validity of our relativistic limits, and for the 
rigor and novelty of the approach, we derive our basic wind 
equations using the covariant formulation of magnetohydrodynamics. 
For the reader unfamiliar with this formalism, we sketch in an 
appendix the derivation of the non-relativistic limits of our 
equations, using the more familiar descriptions of MHD. 
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II. MATHEMATICAL FORMULATION 


a) The Covariant Equations of Motion 

We consider the relativistic magnetohydrodynamic 
equations appropriate for an ionized, highly conducting fluid 
in the presence of electromagnetic fields, following closely 
the formulation of Greenberg (1971) . The space-time metric 
tensor g^^ (Greek indices take the values 0, 1, 2, 3) is 
defined such that 


ds^ = g dx^ dx^ 

with the contravariant four-velocity given by 


(II. 1) 


with U = + 1 (II. 2) 

ds U 

We take the coordinates (x^, x^, x^ , x^) to be (ct, r, 0, <j>) 
in the usual spherical polar sense, and assume that the metric 

A 

tensor g is determined solely by the presence of a spherically 
symmetric body of mass M. If G is the gravitational constant 
and c the speed of light, then we define the "escape velocity" 
c3g such that 


3 


2 

e 


2GM 


rc 


(II.3) 


and we take 


n = 1 - (II. 4) 

We choose the Schwarzschild metrics g^^ = ri , g^^^ = -l/n , 

922 ' ^33 ^ -r^sin^e, and g^^ = 0 for . 



- 6 - 


The antisymmetric electromagnetic field tensor is 

(cf. Landau and Lifshitz, 1971, Chapter 10) 

«3 

^iO ~ ®i ^12 n sin 0 


(II. 5) 


F = + 
^13 


H 2 sin 0 


F 22 = “ r^ sin 6 


where and are space vectors {Latin indices take only 
the values 1, 2, 3). The three-dimensional metric tensor 
is - If is the four-current density, and g is the 

determinant of the tensor then Maxwell's equations have 

the form 


3F 


yv 


3x 


3F 


Ay 


3x' 


3F 


vX 


= 0 


(II-6a) 


dx'- 


/ - g 




3x^ 


- g 


l = -4 


47T ^.v 
_ J 


(II-6b) 


The electromagnetic energy-momentum tensor 
by 



is given 


S 


V 

y 


= — r 

4n L 


yX 4 y aX 



(II. 7) 


We take p to be the isotropic pressure of the fluid, e 
to be the rest energy density of the fluid, and p* to be the 
local rest mass density. The stress-energy tensor T^^ for 
a fully ionized fluid in the presence of an electromagnetic 


field is 
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= (e + p) U'' - pg^'' + S 


P TT^ 


,uv ^ „UV 


(II. 8) 


The equation of motion of the fluid is 


= 0 
?v 


Following Greenberg (1971) , we use the expression for T 
given by equation (II. 8) to write the space components of 

equation (II. 9) as 


(II. 9) 


, » y.A j, y 9 P i T? ^ 

(e + p) U = h ^ - 47 h F 

a 


(II. 10) 


where 


= rrViV _ yV 


h^'' = g 


The equation of continuity of rest-mass takes the form 


(p* u"") = 0 


(II. 11) 


If p* is the local charge density and a the electrical 
conductivity of the fluid, then the invariant form of Ohm's 
law in magnetohydrodynamics is 

= p* - a F^^ U, (11.12) 

We will assume that o is sufficiently large that the electric 
field in the co-moving frame is essentially zero; that is, we 
make the assumption that 

F, =0 (11.13) 

Ay 

With this assumption, the local joule heating is zero, and we 
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may take a polytrope relation between p and p* 

P “ (p*)*^ (11.14) 

The time component of equation (11.9) implies that the energy 
density e is given by [cf. equation (15) of Greenberg] 


e = 


+ 



1 


If we define c3 to be the local speed of sound, then 

s 


e 


2 

S 


gp 


p*c^ 


(11.15) 


and we find that 


& " 

£ + p = p*cMl + - - f - j ) (11.16) 

Equations (II. 1) - (11.16) represent the basic set of fluid 
equations for collisionless plasma. We shall now specialize 
these equations to the problem at hand. 


b) Radial Dependence of Wave Amplitudes in the WKB 
Approximation 

We seek solutions to the above equations which exhibit 
steady-state radial streaming, a radial background magnetic 
field, and time-dependent transverse fluctuations in field and 
velocity which are locally Alfvenic and outwardly-propagating. 
Our procedure is as follows. We decompose our covariant 
equations into radial components and components transverse 
to the radial. The eikonal (or WKB) method (Weinberg 1962) 
is applied to the transverse equations to obtain the local 
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dispersion relation and the wave amplitudes as functions of 
radius, density, radial velocity, and radial magnetic field. 
Once these expressions are known, we can write the radial 
component of the momentum equation and the conserved energy 
flux equation as functions of radius, density, and radial 


velocity alone. 

In this section, we find the local wave amplitudes using 
the short-wavelength approximation, often referred to as the 
WKB approximation. We initially assume that the Alfven waves 
are linearly polarized in the ())-direction, and then generalize 
our results to the case of circular polarization. From Mjzfller 
(1952, Chapter X), we may write the four-velocity of tlje 

plasma as 


= r(i, e, 0, - 


_6B 

sin 6 


•) 


where cB is (dxVdt) and thus represents the contravariant 


component of the radial spatial velocity; c6g is the transverse 


velocity, and 


r = (n - (11.17} 

For transverse Alfvenic perturbations the density and field 
strength perturbations are zero to first order. The conser- 
vation of mass [equation (II. 11) } implies that 

Ip (r"p*re) = 0 (11.18) 

If is the covariant component of the time-independent 
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radial field, then one of the two non-trivial relations con- 
tained in equation (II. 6a) is 

Ip (r^H^) = 0 {II. 19) 

By assumption, the transverse contravariant component H is 
zero. If we take 


* , — w 

r sin 6 


then equation (11.13) implies that 

= 0 

E® = 


where 


6E 


B6H 

n 


<5B H 

r 


( 11 . 20 ) 


( 11 . 21 ) 


From the second non-trivial relation in (II. 6a) we have 


E It ? If (11.22) 

The 0-component of the space part of the equation of motion 
[equation (II. 10) ). is 


_ r-p 2 

47rn - (G + p) = 0 (11.23) 

and the (}>- component of the same equation is 

(e + p) [ r^ i 63 + |I (rrSB)] + r"e 63 

= 4? 1 wit 


(11.24) 
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Equations (11.22) and (11.24) can be solved for the WKB 

amplitudes under the assumptions that; 1) the wavelengths are 

small compared to local scale heights; 2) 60 and 0^ are small 

compared to one, so that we can ignore second order terms in 

60 and 0 . We do not assume that 0 or 0 are small compared 
s © 

to one. Since the solution to these equations is tedious we 
refer the interested reader to Appendix A, and merely quote 
the result. With the above assumptions, our solutions for the 
fluctuating quantities take the form 

60(r,t) = 60(r) exp [i(mt - S(r))] 

6H(r,t) = 6H(r) exp [i(oit - S(r))l (11.25) 

6E(r,t) = 6E(r) exp [i(wt - S(r))l 

where w is angular frequency and dS/dr is the wave number of 
the wave. For notational convenience, we shall not distinguish 
between the full, rapidly varying functions of space and time, 
such as 60(r,t), and their WKB amplitudes, such as 60 (r), which 
are slowly varying functions of space alone. If the meaning is 
not clear from context, we explicitly note the appropriate func- 
tional dependence on space and time. We define the velocity c0 as 

3i 



(11.26) 


In the limit that 0^ is zero, c0^ is simply the local Alfven 
velocity (see, for example, Harris 1957). If we let k = dS/dr, 


then 



the phase velocity of the waves. 


is given by w/k. 
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For outwardly propagating waves, we then find 


6 + 

■ 1 ^ ee.A '“-27) 

G. 

In the limit that 6_ is zero (n = 1) , the phase velocity c3 

e p 

is the proper relativistic sum of the streaming velocity c3 
and the local Alfven velocity cg^. In the absence of streaming 
and in the limit that p* goes to zero, the phase velocity of 
the waves becomes equal to the local speed of light, nc. Note 
that since 0 cannot exceed n [cf. equation (11.17) ] and $ 

A 

cannot exceed one, the phase velocity of the waves cannot in 
any circumstance exceed the local speed of light. The amplitudes 
5H(r), 6E(r), and 66 (r) are simply related by 


SH = 

«E(1 - 6/B„) 

66 H E_ 


(11.28) 


V 


We may verify that the e-component of the momentum equation (equa- 
tion (11.23)] is identically zero for 6E, 6 h, and 6B related in 
this way. The WKB solution for 6E(r) is [see (A. 16) in App. A] 


— g — = const (11.29) 

The corresponding expressions for the radial dependence of 
6H and 60 can be obtained from equation (11.28). In particular, 
we find from equations (11.28) and (11.29) that 

6 ’ r^ 

60^ = const (11.30) 

rMe + T)6^)^ 

a 
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In deriving equation (11,30), we have approximated 1/r^ as 
This is valid in the present context, since we 
are writing an expression for a quantity already assumed to 
be small, and thus need not include corrections to that 
expression which are of order 63^ as compared to the leading 
term. We cannot make this approximation unless our leading 
terms are already small. It is also convenient to express 63^ 
in an equivalent form. We note that the definition of 3^ 
[equation (11.26) ] implies that 



1 - 




47rp*c^ 


(11.31) 


Using the conservation of mass and magnetic flux, we obtain 
from equations (11.30) and (11.31) the form 


63 ^ 


33^(1 - B_") 

const — 

r* (3 + ri3^) ^ 


(11.32) 


In the above discussion, the 4>-polarization was chosen 
because Alfvenic perturbations in this direction can be simul- 
taneously fitted together in a consistent manner over the 
surface of the entire sphere. This is not possible for waves 
polarized in the 6-direction (in particular, note the behavior 
at the poles of a spherical polar coordinate system.). Locally, 
however, the solutions should be valid for arbitrary polariza- 
tions , and in particular in the equatorial plane of our 
coordinate system, both 0- and (|)-polarizations obey the same 
equations. For convenience, in the in what follows we shall take 



the waves as circularly polarized in the equatorial plane, with 
the WKB amplitudes for circular polarization in the same form 
as equations (11.28) through (11.32). The assumption of 
circular polarization has the advantage that for monochromatic 
waves, quantities such as 6E (r ,t) • 6E (r ,t) , 6E(r, t)X5H (r,t) , etc., 
are no longer functions of time. In addition, field strength 
and density perturbations are zero to all orders, with the 
consequence that large-amplitude Alfv^n wave solutions are 
possible (for example, see the treatment by Barnes and Suffolk 
1971) . 

c) The Total Energy Flux 

Having solved for the transverse wave amplitudes, we 

now seek an equation for the conservation of total energy flux 

in the radial direction, including energy flux due to the 

presence of Alfven waves. We first note that the conserved 

mass flux is given bv 
M ^ 

= 4Trr^rp*c6 (11.33) 

The time component of equation (II. 9) is T ~ 0. Using 

the expression for T*^'* given by equation (II. 8), we have 


Ij: {4Trr^(e+p) r^3cn + 6E 6H] } = 0 


(11.34) 


The second terra inside the brackets is just the radial 
component of the Poynting vector. Using equations (11.16) and 
(11.33), we have that the total energy flux is given by 


4Trc p*F3 


-] 


I 


(11.35) 
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This expression includes the energy flux due to the rest mass 
energy associated with the mass flux We consider the non- 

relativistic reduction of equation (11.35) in a subsequent 
paragraph . 

For the moment, assume that we have wind solutions, and 
let cg^ be the radial streaming velocity at infinity, with 
Yoo “ (1 “ Then for a given mass and total energy flux. 


we have 


Yoo = 


E 

C^F 


M 


(11.36) 


The expression in brackets in equation (11.35) is F /c^F„ 

E M 

and is also constant along a streamline. We shall find it 

convenient to write this expression in several ways. First, 

from equation (11,28), 6E 6H = r)6E^/g , with g given by 

P P 

equation (11.27), and from equation (11.29) we have 


6E ^r ^g 
r^g 

ao 

where the subscript "o“ refers to some reference level r^. The 
term in brackets in equation (11.35) becomes 


(1 + 


a - 1 


)rn + 


6E 2 
o 




e (4irp *r 8 c^) e + ne 

ao ^o o o a 


= const (11.37) 


We may also write this in terms of 6g^. Using equations (11,27) 
and (11.28), we have 


6E 6H 


r“ 6g^ H 

r 


(n + gg^) (3 + ng^) 

d d 




a 


2 


2 


(11,38) 
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From equations (11.31) and (11.38), we find that the term in 
brackets in equation (11.35) can be written as 


(1 + _ 3^ )rn + 63 ^ ■ ^ 


= const 


B(1 - B^") 


(11.39) 


®y simple alqebraic manipulation, we can also write equation 
(11.39) as 




(B + nB^)^ 

- B.") ^ 


(11.40) 


63^r 


2 p 3 


(B + nB,) - const 


From equation (11.32) we see that the second term in this 
equation is constant, so that we have for any r the relation 


6 ^ 

^-_^-j) Tri - — (3 + nB^) = const 

8i 


(11.41) 


d) The Radial Equation of Motion 

We have obtained an expression for the conserved 
energy flux in terms of background parameters alone, and thus 
have found implicitly all solutions 3(r) which satisfy that 
equation. To facilitate the imposition of critical point 
requirements, however, we need an explicit differential 
equation for the streaming velocity 3. We may obtain such an 
equation by considering the radial component of equation 
(II. 10). Since we have assumed purely radial expansion with 
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a polytrope relation between p and p*, we may also obtain an 
equation for d0/dr by differentiating any of the above expres- 
sions for the energy flux. We follow the latter course, using 
equation (11.37). The differentiation is a straightforward 
process, although tedious due to the fact that r is a compli- 
cated function of 3 and p* [cf. equations (11.17), (11.26), 

and (11,30) ]. First, we write equation (11.18) for the 
conservation of mass as 


k- ial + 

p* dr r dr 8 dr r 


(11.42) 


In differentiating equation (11.37), we encounter terms 
involving derivatives of p* , and we systematically eliminate 
them in favor of derivatives of d8/dr and dF/dr, using equation 
(11.42). For example, we can easily show from equations (11.19) 
and (11,26) that 


dB, 

dr^ 


- 1 - 


3,^) (- 


dp* 

dr 




which is simply rewritten as 



6 ") 


^8 dr 


+ - 1 ) 
r dr r^ 


(11.43) 


(11.44) 


Proceeding in this fashion, and using simple relationships 
such as dri/dr == (cf. equation (II. 4)], 

differential form of equation (11.37) as 


we write the 
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—inii + 63^) + (nB^^ + 288 ^ + n)r=^68V2] 

= I §|[Tir8s"- r'68MriBa" + 288a - n)/2] (11.45) 

3 ^ 

- ^[Ba'^rd + -~r^) - 2 nr 8 g^ - r^ 68 Mn + 2883 

To obtain a differential equation for 8 alone, we 
differentiate equation (11.17) for T to obtain dr/dr in terms 
of d8/dr. Using equation (11.30) for 68^/ we find after some 
effort that the differential form of equation (11.17) is 


dr 


I? * n/lXe ^ ^ 8 a'>l) 


' fH 4ie‘.°nB - 8a’> ' 6U + 38^^))} (11.46) 

di 

6 2 r ^ 

I {^(1 ^ . 27^1^18(1 - 36 /) . 26//- ne/l . 6 /)]} 

n 3 


3o2 


3 rn2 


We now combine equations (11.45) and (11.46) to write an 
equation for d8/dr alone in the form 


B dr (r,8) 


(11.47) 
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e) The Nbn-relativistic Reduction of the Equations of 
Motion 

To demonstrate the correspondence between equation 
(11.47) and the familiar equation of motion for polytrope 
winds, we consider this equation in the limit that 3 and 
are small compared to one. We do not assume that 3^ is small ^ 
compared to one. Neglecting third order terms and higher in 
small quantities in equations (11.45) and (11.46) [e.g., 363^? 

3^^63^, etc.], we find that in this limit d3/dr is given by 


r 

3 dr 




63^ 

“ 213+0 

Cl 

, 53 ^ 

41071 

a 


(1 + 3 ^^ 


(1 - 3 ^^ 


) (3^ + 33) 


) (3^ + 33) 

Cl 


(11.48) 


We have kept terms of the form 3 63V (3 + 3^) (which at 
first appear to be third order) in equation (11.48) because of 
the possibility that 3^ << 3» Equation (11.48) can also be 
derived by more familiar techniques. For the convenience of 
the reader who is not at ease with the mathematical formalism 
used above, we derive equation (11.48) in Appendix B using the 
standard MHD equations when not only 63 and 3g but also 3 and 
3^ are small compared to one, with 3^ unrestricted. In the 
limit that 3 is small compared to one, equation (11.48) 

3i 

iTsducGs to th 0 oqiistions of motion ussd by Aisziraki 3.nd. Coutuiriox* 
(1971) and Belcher (1971). If 63 is zero, we obtain the 
standard form for the equations of motion of polytrope stellar 


winds o 
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We also demonstrate the correspondence between our 
expression for the total energy flux F„ and its more familiar 
forms. As above, we expand F_ assuming that B and B are small 
compared to one, and neglect third order terms in small 
quantities. In this limit, from equations (11.35) and (11.39) 
we have 

Ffi = c^Fj^ + 47rr^{cB[j c^p(B^ + 6B^) + - ^ p 

(11.49) 

In the limit that B is a small compared to one, this 
expression for reduces to that given by Belcher (1971, 
equation 26b) except for the rest mass energy flux term, c^F^j. 

For future convenience, we rearrange terms in equation (11.49) 
to obtain the form 


3 + 3 

- j c^pB ^ ] + pc*6B^ —} 

1 - 3,^ 


Fv. = I f t e." C 


2 *"e 


B^ + 2B(1 - 6,V2) 

+ 6B^ ] 

3(1 - B,=") 


(11.50) 


where C is given by 




(11.51) 


Note that the expression in brackets in equation (11.50) is 
Yoo “ ^ situation that B and B^ are small compared to one. 
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III. NUMERICAL SOLUTIONS TO THE WIND EQUATIONS 


In the preceding section, we have set up and formally 

solved the equations of motion for the dynamic expansion of 

a stellar corona. We now consider solutions to these equations 

which satisfy given initial conditions at some reference level 

r close to a coronal base. In this section we sketch a 
o 

general numerical algorithm for obtaining wind solutions for 
given initial values. We exhibit a limited number of these 
full numerical solutions, and note some of their characteristic 
features. In Section IV, this qualitative information will 
enable us to obtain approximate analytic solutions to the 
critical point equations over some ranges of initial values. 

We choose a value for the poly trope index a, and specify 

at some reference level r^ the escape velocity 

the Alfven velocity the sound velocity and the 

velocity perturbation c60^. If we also choose a value for 

eg , the radial velocity at r , then this value of g along 

with the set of initial conditions in the combination 

{g , g , 6g ^/g g ^/g is sufficient to determine 

eo' ao o eo so eo 

g as a function of r. This function is not necessarily a wind 
solution to equa:tion (III.l), as we have imposed no critical 
point requirements. To determine g (r) . we first note that 
equation (11.39) for the conserved energy flux divided by the 

mass flux is 


3 ^ 

(1 + — s ) m + 6g^r^ 

' a - 1 


(n + gg^) (g + ng_) 

ci Q 

g(l - g^^) 

d 






(III.l) 
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where the subscript ”o" refers to the reference level. All 
quantities on the right hand side of equation (III.l) can be 
computed using our initial parameter set and 3^* On the left 
hand side, for a given r, we easily compute the local escape 
velocity using equation (II. 3) in the form 3 ^ where 

we define Z as 

z = |- (III. 2) 

o 

Of course, ri at r is then 1 - 3^^. We now guess a value 3 
for the radial velocity at r, and check to see if this value 
satisfies equation (III.l). This process is complicated by 
the fact that r at r depends on 63^ at r, which in turn depends 
on r at r in a complex way [cf. equations (11.17) and (11.30)]. 

To obtain an initial estimate of T we approximate F at r by 

(n - 3^/n) Using the conservation of mass [equation (11.18)], 

we estimate p*/p * at r to be F 3_/F3Z^. Given an estimate of 

O O G 

P*/Pq* at r, we obtain an estimate of 3^ at r using equation 
(11.26) in the form 

1-3 ^ 

3_ = [1 + — ^ Z**]"^ (III. 3) 

® 3 " *'o 

ao 

We then compute a first estimate of 63^ at r using equation 
(11.32). We now improve our estimate of F at r by using this 
estimate of 63^ in equation (11.17). We then compute a new 
estimate for p*/p^* using the conservation of mass, and a new 
estimate for 3^ and 63^ using equations (III. 3) and (11.32). 
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This leads to a better estimate for F, and so on. This 

I ^ 

process is iterated until the n estimate of T differs from 
the (n-1)^^ estimate by less than one part in 10^®^ 

square of the sound velocity at r [cf. equation (11.15)] is 

then 6^-6 ^(p*/p Given our guess for 3 at r, we 

s so o 

have thus computed values for n, 63^/3^/ T, and 3^^ at r. 

We now check our guess for 3 by using these quantities to 
determine if equation (III.l) is satisfied. If not, we keep 
guessing until we find a value of 3 that does satisfy 
equation (III.l). Thus, for a given 3^ and our initial values, 
we can determine 3 as a function of r . 

We find the wind solution (if it exists) by choosing 3^ 

such that 3 as a function of r passes through the critical 

point of the differential equation (11.47). For a fixed 

initial value set and variable 3^/ the critical point 3^) 

is determined by the requirement that , 3^) and 

c c o 

^'(r ,3,3) simultaneously vanish (we have explicitly noted 
c c o 

the dependence of ^and^ on 3^) . In addition, (I’^f 3^) must 
also lie on the solution 3 (r) which satisfies equation (III.l)? 
this requirement imposes a third condition of the form 
^{r^, 3^, 3^) =0. We thus have three transcendental equations 
which determine the three quantities r^, 3^> and 3^^, and thus 
the wind solution 3(i^) • We refer to the combination 
(Zc, 3^/3g^, 3 q/3^q) / where solution set for 

the initial value set (3^^, 3^^, ^so'/^eo'^* 
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determine the solution set for a given initial value set in a 

manner similar to the procedure used to determine 3 (r) from 

equation (III.l). That is, we guess a solution set 

(Z , 3 /3 , 3/3 ) compute the quantities T , p */p *, 3 , 

3 ,63 and n exactly as described above (the subscript "c” 

SO o o 

refers to the critical point) , and check our guess by seeing 

if y, and simultaneously vanish. Thus, given our initial 

value set and the additional critical point requirement, we 

determine the solution set and the wind solution for 3 as a 

function of r. For < r < r„, we look for solutions to 

o c 

equation (III.l) in the range 3 < 3 < 3 / and for r > r , we 

o c c 

look for solutions in the range 3 > 3 • 

In practice, of course, the transcendental equations which 
determine a solution set for a given initial value set are 
involved, and it is impossible to guess a correct solution set 
for these equations £ priori . Computationally, our procedure 
for finding solution sets always begins with a known solution 
for a given initial value set. To find a new solution set for 
a different initial value set, we slowly vary one of the initial 
values and numerically follow the solution set into new regions 
of solution space, starting from the known solution set. Since 
our known solution set originally began with the Parker solutions 
to the stellar wind equations, all of our numerical solutions 
are in a real sense analytic continuations of Parker solutions 
into new regions of initial value space. 

Using this algorithm, we have computed full numerical 
solutions for a large range of initial values, and we display 
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soitie of these solutions to indicate their general characteristics. 

We take to be the quantity C as defined in equation (11.51) 

evaluated at the reference level. We recall that for Parker 

polytrope winds to exist, we must satisfy the inequality 

^ ft 2 /g ^ >1, Since we are considering only conditions 
a-1 ^so ' ^eo 

under which thermally driven winds do not exist, the quantity 

is always positive. To minimize thermal effects, we take the 

polytrope index a to be 5/3, the adiabatic value. We plot our 

solution sets as functions of with the remaining 

initial values held fixed for a given curve. Figure 1 gives 

values of Z^, ^o'^^eo functions of 66^ /^eo 

for a value of log fixed at - 4.0 (c3 = 30.0 km/sec} and 

a value of log equal to - 2.7 = 598.6 km/sec). The 

curves labeled A, B, and C are for values of log (3gQ /^eo ^ 

equal to - 1.0, - 1.4, and - 3.0 (corresponding to sound 

velocities of 189.3 km/sec, 119,4 km/sec, and 18.9 km/sec, 

respectively) . We can determine ^ function of 

66 ^/6 and this is also plotted in Figure 1. 

^o ' eo 

For these same initial values, we plot in Figure 2 the 
constants of motion F„ and - c"F„, and the quantity F^/o"f„- 1, 
which is Y„“ considering the scales for the mass and 

energy fluxes, we must remember that we have not yet completely 
specified conditions at the reference level. Up to this point, 
we have used the reference density and field only in the com- 
bination and the reference level and stellar mass 

only in the combination M/r^. These combinations are sufficient 
for the solution of the critical point equations, but to express 
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our mass and energy fluxes in grams/sec and ergs/sec, we must 

in addition specify M (or r^) and p^* (or H^) . We choose to 

consider our remaining independent variables as the stellar 

mass M and the n\imber density N (where N is the rest mass 

density divided by the proton rest mass) . Once these variables 

are specified, we have that r^ in centimeters is given by 

3 X 10® if M is measured in solar masses, and H ^ in 

eo ' o 

gauss squared is given by 0.019 N 6 /(I - 3 if M is 

o ao ao o 

measured in number per cubic centimeter. It is clear from 
equations (11.33) and (11.35) that both and are propor- 
tional to To obtain absolute units in Figure 2, we choose 

M to be one solar mass and to be 10®cm“®. Mass and energy 
fluxes for other choices of M and can be obtained by scaling 
according to For our present choices, with log 6^^ = - 4.0 

and log 3__ = - 2.7, we have r^ = 7.5 x 10^® cm and H =0-14 
gauss. For reference, we note that for the solar wind, mass 
and energy fluxes are on the order of 10 grams/sec and 
3 X 10^ ’ergs/sec, respectively. 

There are characteristic features of the curves presented 
in Figures 1 and 2 that are common to all of the numerical 
solutions obtained. First, the curves are not too sensitive 
to the temperature parameter , particularly when this 

ratio is small compared to one. Note that for moderate wave 
amplitudes, the energy flux F^ - is essentially independent 

of temperature. Second, there is a broad range over which our 
solutions exhibit power law dependencies on the wave amplitude 
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66^/3^^, In this range, 3 and are essentially constant, 

o eo ^ c eo c 

F„ and 3 /3 go as 63 '* /B i F-c. “ goes as and 

M o' eo ^ o' eo E M ^ o eo 

- 1 goes as ^ ^ future reference, we call 

this power law regime the intermediate wave amplitude range . 

As we move toward lower wave amplitudes, we eventually encounter 

an abrupt, low-amplitude cutoff in our dynamic solutions, at 

which 3 q and the mass flux abruptly fall to zero and the energy 

flux is well behaved, so that Y«. “ 1 goes to infinity. This 

abrupt cutoff occurs as approaches 43 5 V27 from 

o eo cto o 

above (for the curves A, B, and C in Figure 1, the cutoff occurs 
at = 5,1 X 1 Q-®, 1,31 X 10"^and 1.47 x iQ-s, 

respectively) . For wave amplitudes below this cutoff, wind 
solutions do not exist, but static solutions are possible, as 
we shall see later on. As we move toward higher wave amplitudes, 
we encounter two different phenomena, depending on whether 3 

ao 

is substantially greater than or smaller than 3 . In the case 

eo 

8^0 ^ ^eo shown) , the power law behavior of the intermediate 

range holds until becomes comparable to one, at which 

point B^/BgQ and approach one. If we move substantially 

beyond this limit, r^ becomes less than r and 3 exceeds B 

c o o eo 

In the case 3^^^ < 8 ^^, the power law behavior holds until 

‘^^o^'^^eo^ approaches As is evident from Figure 2, 

Fg - c^Fj^ no longer increases as 83 ^^ 73 ^^^ beyond this point, 
and in fact begins to decrease, with r^/r^ increasing. Thus, 
as ^B^^/BgQ^ becomes comparable to the smaller of 1 and 
3ao/8gQ t our simple power law dependencies disappear, and we 



- 28 - 


refer to this regime as the strong wave amplitude region . 

Having determined the value of 6 /8 for a given set of 

o eo ^ 

initial values, we can compute 8(r), p*(r), and so on, using the 
algorithm previously described. In Figure 3 , we show typical 
wind profiles for solution points on curve A of in Figure 1, 
as marked. The profiles 1 through 5 in Figure 3 are for 

^so"/^eo" " ~ 

and - 5.32, respectively (corresponding to c6B^ = 75.4, 23.8, 7 . 5 , 

2.7 and 1.3 km/sec, respectively). The last amplitude in this 

series is below the cutoff for dynamic solutions, and the 

corresponding profiles represent static solutions, as discussed 

in Section V. Figure 3a gives profiles of p*/p * and B/B as 

o eo 

functions of r/r , and Figure 3b shows g and 6 B/B as functions 

of r/r^. We note that <SH/H is equal to - 6 B/B if B « 1, so 

that the profile of BB/B^ is essentially also that of 6H/H. 

Figure 3c gives the wave amplitude 6B divided by B and the 

eo 

ratio of the transverse velocity to the radial velocity, 6B/B. 

In Figure 3d we plot the Alfvenic Mach number B/B , and the 

total wave energy flux A, normalized to its value at r . The 

o 

energy flux A includes all terms in the total energy flux which 
are proportional to the squares of wave amplitudes. The vertical 
line on each curve marks the location of the critical points. 

Except for the strong wave amplitude case {curve 1) , the 
radial profiles exhibit rapid decreases in density at a radial 
distance which in the absence of waves would be the top of a 
static atmosphere, with correspondingly rapid increases in B , 

Si 



3, and 53* Beyond this point, the density begins to fall off 

as 1/r^ , the velocity increases slowly, approaching the constant 

c3 and 63 ^n.d 3 eventually begin to decrease, approaching 

zero far from the star. For future reference, we point out some 

of the prominent characteristics of the profiles corresponding 

to the intermediate range of wave amplitudes (curves 2, 3, 

and 4). First, because of the rapid decrease in density, 3^ 

at the critical point is small compared to 3 and 63 there, 

even if 3 were much greater than these quantities at the 
S 

reference level. Second, from Figure 3d we see that 3 is 

much less than 3 both at the reference level and the critical 

a 

point. Finally, and most importantly, we consistently find 

that 63 and 3 are approximately equal at the critical point, 

with a value close to the escape velocity there, regardless 

of their initial values at the reference level. Because of 

this property, the different profiles of 63 and 3 in Figures 3a 

and 3c are hard to separate; in general, the lower the values of 

6 3 and 3 vesry close to the star, the higher their values very 

far from the star. From Figure 3d, we see that the transfer 

of wave energy flux to kinetic streaming energy flux takes place 

gradually, with no steep gradients, and with essentially all of 

the wave flux at r^ going to streaming energy far from the star. 

As 66 ^/3 ^ decreases toward the low amplitude cutoff, 

o eo 

the density profiles become increasingly rarefied, with the 
Alfv4n velocity at the critical point approaching c. Beyond 
the critical point, 63 and 3^ require much greater distances 
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for significant decrease, and B approaches its limiting velocity 

much more slowly. As approaches the smaller of 1 and 

^ao^^eo ' density profiles fall off less rapidly, and the 

gradients in 3 and 6^ are correspondingly diminished. Above 

this point, either Q becomes comparable to 6 (if 6 > 6 ) 

o eo ' "^ao ^eo 

or 3 becomes comparable to 3 (if 3 <3 ) . 

^ 2 lc ao eo 
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IV. APPROXIMATE SOLUTIONS TO THE CRITICAL POINT 
EQUATIONS IN SPECIAL CASES 

a) Basic Assumptions 

With experience based on the full numerical solutions 

presented above, we seek approximate analytic expressions for 

our solutions (Z , 3 /B ) over various range of 

initial conditions. At the outset, we limit ourselves here 

to cases for which 6 << 1 and 3 << 1. Additional limita- 

eo ao 

tions on the ranges of our initial values will appear as we 
proceed. To obtain analytic solutions, we make reasonable 
assumptions as to conditions at the critical level and 
reference level, and we list the assumptions below. After 
finding solutions, we a posteriori check the validity of these 
assumptions. 

Assumption As 3 << at the critical level and at the 
reference level. 

With this assumption, 3 must be small compared to one at 
the critical point. Since we have already assumed that 3^ is 
small compared to one, the non-relativistic equation of motion 
(11.48) is appropriate, and our critical point occurs when both 
the numerator and denominator of this equation are zero. 

Dropping terms such as 63^3/3^ (third order small by assumption) , 

a 

the critical point equations become 


ec 


S3. 


- 3 


sc 


(1 + 3 _") 


ac 


4 


2 


4 


(IV. 1) 
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and 


63 ^ 

B 2 - 3 2 (1-3 

c sc 4 ac 


(IV. 2) 


Note that we have not assumed 3 to be small compared to one. 
We may combine equations (IV.l) and (IV. 2) to obtain 



5 ^ (1 - 3 23 2 

ec ac . ^ac 

A * 


<1 * 63^^) 


1 + s 


sc 


ac 


(IV. 3) 


For the integrated form of the equation of motion, we use 
equation (11.41). To second order in small quantities, we 
have 


3 ^ + 63 ^-3 ' 

*^ec 


2 

a-1 


3 " = 

sc 


eo 




(IV. 4) 


To find 63^ at the critical point, we use equation (11.30). 

Under our present assumptions, neglecting terms of order 

63_^B /3__, etc. , we obtain 
o O 30 


63 


c 


2 



ac 



(IV. 5) 


Assumption B; 3^^ << [3^^^ + - 1] 

With this assumption, we may drop 3^^ on the right-hand 
side of equation (IV. 4). We take to be the quantity ^ as 
defined in equation (11.51), evaluated at the reference level. 
As we have noted previously, we only consider situations in 
which > 0. Then equation (IV. 4) becomes 
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ec 


- 


a - 1 


sc 


= e 


eo 


(IV. 6) 


Assumption C: << 3^ at the critical point. 

This condition is one of the novel properties of wave 
driven winds. Xn thermally driven winds, of course, is on 

the order of at the critical point. The dynamical rationale 
for this condition will become clear when we consider the 
transverse velocity c6B at the critical point. With this 
assumption, we may drop the terms containing B^^^ from the 
right-hand side of equation (IV. 3) and from the left-hand side 
of equations (IV. 1) and (IV, 6). As we shall see below, the 
only time this is not justified is when B^^ is extremely close 
to 1. Equation (IV. 3) is thus taken to be 


B 


2 

C 





(1 . 


(IV. 7) 


Equation (IV. 1) is 


6B 


2 

c 


ec 


1 + B 


ac 


(IV. 8) 


and equation (IV. 6) becomes 


^c' ^^c' - ^c' = ^eo^ ^o 


2 


2 


2 


2 


(IV, 9) 
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Using equations (IV. 7) and (IV. 8) in (IV. 9), and remembering 

that we obtain an equation for Z as a function 

of B and C 
^ac 

7 + 58 ^ , 

Z = — — — (IV. 10) 

4(1 + 8 ? 

"^ac ^o 


We now solve for 8 in terms of the conditions at the 

ac 

reference level by equating the expression for 68^ from 
equation (IV. 5) with that in equation (IV. 8). If we define 
the variable Xq 


X = ^ 

Xq 4 


^ao^o 




eo 


(IV. 11) 


where is defined by equation (11.51), then the resulting 
equation may be written as 


X, 


( 12 ) 3 (1 + 

ac 


= 1 


(IV. 12) 


Equation (IV. 12) determines 8 as a function of the initial 

ac 

conditions. Given 8 , we immediately have Z from equation 

a c c 

(IV.IO) , and from (IV. 7) we find that 


B 

B 


c 

eo 



- 6ac » 


(1 + B ) 

ac 


H 

h 


(IV. 13) 
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To complete our solution, we need an expression for Pq/^qq’ 
From the conservation of mass, we have 






(IV, 14) 


Using equation (III, 3), we obtain as a function of 3^^, 

and Z^, so that equation (IV. 14) becomes 



(1 - 


ao 


(1 t 3^^^) 


2Z_Vz 


ac 


(IV, 15) 


At this point, we have completely solved the critical 

point equations under the limiting assumptions A, B, and C. 

Given the initial values (3^^, 0^^, ^so"/^eo"^ ' 

we compute the parameter Xq using the definitions in equations 

(11.53) and (IV. 11). Given x^, we then find that value of 3^^ 

o 

between zero and one which satisfies equation (IV, 12) [if it 

exists]. The solution for (Z^, 3^/3^^, ^o'^^eo^ follows 

immediately from equations (IV. 10), (IV. 13), and (IV. l5). To 

demonstrate the range over which these expressions are valid, 

we plot in Figure 1 the analytic solutions along with the full 

numerical solutions for 3 /3 and r /r as functions of 6 3 

with 3 ^/3 ^ = .1. The analytic results from the above equations 

so eo 

are indicated by dots. Our approximate analytic results are indis- 
tinguishable from the numerical results at low and moderate wave 
amplitudes, but at high wave amplitudes there is considerable dis- 
agreement between the two. In the following sections, we consider 
these regimes in more detail, in particular the low-amplitude 
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cutoff in our expressions, and the differentiation between 
"moderate" and “high" wave amplitudes. 


b) The Low-Amplitude Cutoff 

We examine the regions in initial parameter space near 

the abrupt, low-amplitude cutoff in 6 and $ when considered 

o c 

as functions of (e.g., the left-most segments of 

the curves in Figure 1 and 2). From equations (IV. 7) and (IV. 15), 

it is obvious that these cutoffs occur as 8 approaches one 

(or, equivalently, as approaches zero) . From equation 

(IV. 12) , 8 becomes one when is one. If we take 1-8 ^ 

du o '■^ac 

to be small and expand equation (IV. 12) to first order in this 
small quantity, we obtain 


1 - 




(IV. 16) 


Using this form in equations (IV. 10), (IV. 13), and (IV. 15), 
we have 


Z 

c 



8 


eo 




(IV. 17) 






1 ) 




As we show below, these solutions are valid in the ranges of 
initial parameters satisfying 

6 


{— ) 8 


1 2 

<< Y 


4 


1 << 1 
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To find the constants of motion for these solutions, we 
evaluate the expressions for F.j^ and [equations (11.33) and 
(11.35)] at the reference level. Since 3^ and 8^^ are much 
less than one at r^, we use the non-relatxvistic form for F^ 
given in equation (11.50). If we insert the solutions (IV. 17) 
into (11,50), we find that the third term in (11.50) is much 
larger than the first two. The energy flux F^, is to an 
excellent approximation given by 


“Vl'o* (IV. 19) 

The second term in equation (IV. 19) is just the radial 
component of the Poynting flux associated with the Alfven 
waves at r^. Note that all of the wave energy flux at r^ 
appears at infinity in the form of streaming energy flux. The 
mass flux is simply 4wr^^ p^*c3^, with 3^ given by equation 
(IV. 17). The quantity Y„“ ^ given by the third term 

inside the brackets in equation (11.50) , and since 3^ << 
we have 

Y - 1 = 63 (IV. 20) 

'00 

Evaluating this expression for the solutions given by equation 
(IV. 17), and remembering that these expressions are valid only 
when X very close to one, we have to a good approximation 


Y 


00 


V2 , 

«eo^o 



1 


(IV. 21) 
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The restrictions on the range of applicability for the 
solutions in equations (IV. 17) arise from the requirement 
that our solutions must satisfy a posteriori the assumptions 
A, B, and C under which they were derived. The requirement 
that I Xq“ 1 I << 1 easily satisfies both A and B. Assumption C 
is more difficult to satisfy, for the following reason. Using 
(III. 3) and assuming (IV. 16) and (IV. 17) to be valid, we find 
that the density at the critical point is 

^ = 2(X„ - 1) '^''•22) 

Since that goes to zero as 

S vJ SO O O S Cr 

(Xq - 1)*^”^ as approaches one. However, if a is less than 2, 

according to equation (IV. 17), 6 ^ goes to zero faster than this, 

as (x - 1). As a consequence, even though g ^ is usually much 

smaller than in the limit of extremely low density p^* (x^ 

extremely close to one), 3^^^ becomes comparable to 3^^* Taking a 

to be 5/3, we can derive the lower limit in inequality (IV. 18) from 

equations (IV. 17) and (IV.22) , and the requirement that (x^-l) 

be large enough to insure that $ ^<<8 For all cases 

s o o 

considered here, the lower limit given by inequality (IV. 18) is 
extremely small, so that Xq must be very close to unity before 
assumption C is violated. For example, none of our numerical 
solutions near cutoff even approach the lower limit in inequality 
(IV. 18), so that the region below this limit and above Xq ~ 1 


is very narrow. 
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For completeness, however, we give (but do not derive) 
the solutions appropriate to the situation when the lower limit 
in inequality (IV. 18) is badly violated (e.g., >> 3^^)« I 

this circumstance, we assume that is so close to one that 

the second term on the right hand side of equation (IV. 3) 
dominates our expression for • Iri this approximation, if 
we keep terras in equation (IV. 6) only to order ^ 

2 y 

particular, neglecting corapared to ^^”^ac ^ 

find that 
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3 
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eo 


o 


eo 
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Vz 


eo 

3 a 
so 
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1 ) 


This solution is valid in the limited region of initial 
parameter space defined by 


(IV. 23) 


« < X, - 1 « 


3Vi 


(IV. 24) 


The upper limit in inequality (IV. 24) derives from the require 
ment that the second term on the right hand side of equation 
(IV. 3) be much larger than the first. The expression for 
in this range is still given by equation (IV. 19). The mass 
flux is proportional to 3^, as before, and ^ 
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- 1 = 


3g 3 

o ao so 


27 e 


eo 


<Xo-l) 


(IV. 25) 


In deriving the lower limit in expression (IV. 18) and the 
equations (IV. 23), (IV.24) , and (IV. 25), we have for con- 
venience assumed a = 5/3. All other equations in this paper 
are for arbitrary values of a. 

The above expressions [equation (IV. 17) in the range 
(IV. 18), and equation (IV. 23) in the range (IV.24)] exhibit 
all of the cutoff properties of the numerical solutions in 
Figures 1 and 2. The mass flux goes to zero as ~ 1 

goes to infinity as Xq"^!/ and the energy flux is well-behaved 
and approaches a well defined limit at cutoff. If we take 
Xq= 1 , we see from equation (IV. 19) that the minimum energy 
flux we can get from the system in the form of a wind is 


min 


8 

27 


^eo^o 


[ 47rr ^c8 


eo 


H ■ 
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8tt 


(IV. 26) 


The term in brackets is the energy flux that would arise if 
the reference magnetic field energy density were convected at 
the escape velocity cB^^. We let denote this energy flux; 


= 


4TTr 

o eo 


H 2 
Q 

8tt 


(IV. 27) 
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c) The Intermediate Amplitude Case 

We now consider regimes where our initial wave 
amplitudes are well above the cutoff point, e.g., when the 
parameter Xq defined in equation (IV. 11) is well above one. 
In this regime, 6 as given by equation (IV. 12) is small 
compared to one, and we have 


ac 


(i)* 


Equations (IV. 10), (IV. 13), and (IV. 15) then give 


4?, 


eo 


£ »/2 
^o 
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1 27 ^ 7 ® 4 54 -V. 
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^eo 


(IV. 28) 


(IV, 29) 


The numerical factor in the expression for is approxi- 

mately 3.55. These expressions are valid in the regions de- 
fined by 


4 63 ^ 3 

I 7 3^^ ^ < MINd, ^) (IV. 30) 

27 ao o 2 ^ 

eo 

The symbol MIN(a,b) denotes the smaller of the quantities 

a and b. The lower limit here insures that we are well 

away from cutoff, with y >>1 and 3 <<1. The upper limits 

o ac 

stem from assumption A (3 <<3 ) and from assumption B 
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* Assumption C is easily satisfied, and in 
particular 



14.8 



so 


6 ^ 

eo 


6B„ «/3 

( 9 .) 

eo 


(IV. 31) 


We point out that if ‘SBq/B^q<< 1 and 6BQ/BgQ<< 1, then the 
upper limit in inequality (IV. 30) is automatically satisfied. 

We again evaluate the constants of motion at the reference 
level. Within the range defined by expression (IV. 30), the third 
term is equation (11.50) for F^, at r^ is again much larger 
than the first two. For the intermediate case, our expression 
for Fg is simply the sum of c^Fj^ and the wave energy flux at r^ 

[cf . equation (IV. 19) ] . Again we note that all of the wave energy 
flux at r^ appears at infinity in the form of streaming energy 
flux. The mass flux F^ is 47rr^^ p^*cB^, as before, and 

3 3 

- 1 = 0.28 B._ f (IV. 32) 

6B ^ ao o 
o 

To obtain some idea as to the allowed range of the constants 
of motion, we let vary between the limits imposed by 

inequality (IV. 30). Within these limits, F^ - c^Fj^ varies 
over the range 


8 

27 


eo 
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2W^ MIN(1, 
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ao 


(IV. 33) 
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and is proportional to mass flux varies 

between 


V2 '^o 

^o ;T < < 


W B ^ 

° MIN(1, 


g 2^ 7/2 £,2 

eo ^o 


e 


ao 


(IV. 34) 


and is proportional to *5^Q'*/^eo'* ' term W^/c^ [cf. equation 

(IV. 27)] in expression (IV. 34) is the mass flux which would arise 
if the equivalent mass density H^^/8 ttc^ were convected outward 
from r^ with velocity ^8^^. The quantity Yoo “ 1 varies between 
the approximate limits 


2.0 8^0^0^" > Yoo - 1 MAX(1, •^) (IV. 35) 

6 o 

The symbol MAX(a,b) denotes the larger of the quantities a 
and b. The upper limit here has been assumed to be much less 

than one, so that Yoo ” ^ ~ *^^00 ^eo 

c6 < then inequality (IV. 35) can be written as 
eo 


2.0 



/cV 


eo 


> V > 
00 



V MAX(1, 
eo 



(IV. 36) 


In the intermediate range of wave amplitudes, we thus expect 
to find streaming velocities at infinity in the range given 
by inequality (IV. 36), with an inverse dependence on 
All of these functional dependencies on <^Bq/Bqq agree with 
the power law behavior we expect on the basis of the full 
numerical solutions for the intermediate range of wave 
amplitudes (Figures 1 and 2) . 
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d) The Strong Amplitude Case 

In the strong amplitude case, when our initial wave 
amplitudes are greater than the upper limit given by inequality 
(IV. 30), we have not been able to derive explicit analytic ex- 
pressions. In the case that > ^eo' from equation 

(IV. 29) that 3^ becomes comparable to 3^^ as the wave amplitude 
approaches the upper end of the intermediate amplitude range. 

This is an unrealistic situation for a reference level near the 

coronal base. In the case that 3 < r the wind becomes super- 

ao eo 

Alfvenic at the critical point as exceeds ^^o'^^eo' 

Assumption A above is no longer justified, and our intermediate 
range solutions are inappropriate. In the strong amplitude 
regime, we note that in the expression for the energy flux given 
by equation (11.50), we can no longer neglect the - 1/2 
term, since even at the upper end of the intermediate amplitude 
range it is becoming comparable to the third term. The quali- 
tative behavior to be expected would be a decrease in Fj.-c^Fj^ 
as <S3_^/3^rt^ increases beyond 3 /3 / since we are now 

subtracting two terms of comparable magnitude, and indeed this 

is the behavior our numerical solutions for 3 exhibit 

ao eo 

in the strong amplitude domain (cf . Figure 2) . Other than 
qualitative statements, however, we must rely on our numerical 
solutions in this range of wave amplitudes. 
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e) The Alfvenic Critical Point 

In the niameric and analytic solutions obtained above, 
the energy flux associated with the waves decreases as we move 
away from the coronal base, with a corresponding increase in 
the streaming energy flux associated with the wind (cf. 

Figure 3d) . A physical scale of obvious importance in this 
process is the radial distance from the star at which a signif- 
icant fraction of the initial wave energy has been transferred 
to streaming energy. It is easily shown that this distance 
scale is on the order of the Alfvenic critical distance r^ (by 
definition, r is that point at which the radial streaming 

Si 

velocity is egual to the local Alfven velocity) • For example, 

consider equation (11.37) for the conserved energy flux divided 

by the mass flux times c^. The second term in this equation 

is the wave Poynting flux divided by the mass flux times c^ . 

As lonq as & <<3 (i.e., the streaming is sub-Alfvenic at r ), 

this term will decrease by about a factor of 1/2 in going from 

r to r . Of course, there are other wave terms in equation 
o a 

(11.37) in addition to the Poynting flux term, but similar con- 
si dera'*"- ions apply as to the scale height over which these terms 
show significant decrease. We see from our numerical solu- 
tions (e.g.. Figure 3d) that 2^ (i.e., increases as 

6B /B decreases. In the following paragraph, we derive 
o eo 

approximate expressions for appropriate to the intermediate 
range of wave amplitudes. We have not been able to derive such 
expressions for the strong amplitude and cutoff regimes, but 
from our numerical solutions, it appears that is close to 
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one for strong wave amplitudes (if 3 < B ) , and approches 

d.O 30 

infinity as <58^/8^^ approaches the low amplitude cutoff. This 
circumstance again casts doubt on the physical significance 
of solutions just above cutoff. 

We now seek an approximate expression for Z in the inter- 
mediate range of initial wave amplitudes. We recall that in 
this range, 8 and 8_ are always small compared to one, with 

A 

8q << From equation (11.30), we thus have, for any Z, 

3 ^ 

«6“ = 3^ — (IV. 37) 

° ®ao (6 + 

a 


and from equation (III, 3) and the conservation of mass 


8 


a 


^ao . 8 > Vi 
Z ^8^> 


(IV. 38) 


We also have 8^ << 8^^, so that equation (11.41) is, for 
any Z , 


•5 (8 + 8 ) 
8 ^^ - 6 ^ - 68 ^ - .. 8 „^ + 268 ^ - ® 


a - 1 ^s 


8 . 


= 8 C 

eo ^o 


(IV. 39) 


We now evaluate equations (IV.37) through (IV. 39) at Z . From 

A 

equation (IV. 38), we have that 8 at Z is given by 


8 


a 



(IV. 40) 
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Using this expression in (IV. 37), we obtain for at Z 

a 

63 2 p 

63^ = — ^ (IV. 41) 

From equation (IV. 20), we see that 63^ at Z is just 3 ^/8. 

Using these expressions in (IV. 39) and neglecting the sound 
velocity c3_, we obtain a fourth order equation for Z of 

S 3 

the form 
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eo 


2 


Z 


a 




(IV. 42) 


For a rough estimate of Z , we note that expression (IV. 36) 

3 

indicates that 3^^ is significantly greater than except 

at the upper end of the intermediate amplitude range. If we 
neglect terms involving 3^^^ in equation (IV. 42), we find that 

3 ^ 

z = _a^ — 8_ (IV. 43) 

^ 3 ^ 3.3 ^ 

o 


From equation (IV. 40) and (IV.43) , we note that 3^ at Z is 
33^^/8. Using equation (IV. 29) for 3^ and (IV. 20) for 3o^,, 
equation (IV.43) becomes 

Z = 0.80(^^ _®2_) (IV. 4- 

^ ^eo 63 " 

o 

Assuming that 3 < 3 , we find that in the intermediate 

’ ao eo 

range of wave amplitudes defined by expression (IV. 30), Z^ varies 




The lower bound here is suspect, as it occurs at the upper 
end of the intermediate amplitude range, and in this region 
we are not justified in neglecting the terms in equation 

(IV. 42). As long as we stay away from the upper limit in in- 
equality (IV. 30), however, equation (IV. 44) for is reasonably 
accurate, as may be verified by comparison with the full 
numerical results in Figure 3d. 
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V, STATIC ATMOSPHERES 


We have found in previous sections that for initial 
wave amplitudes below a certain minimum value defined by in- 
equality (IV. 30), wind solutions to the equations of motion 
no longer exist. Since the mass flux of the dynamic solutions 
decreases to zero as the initial wave amplitude decreases 
toward the minimum value, it seems plausible to expect that 
below this minimum value only static atmospheres occur. This 
is indeed the case, as we shall demonstrate. 

Even though the development in Section lie above assumes 
throughout that a dynamic expansion exists, we may extract the 
equations appropriate to the static situation by taking the 
limit 3 0. We do not use equation (11.39) for this purpose 

because of the singularity as 3 -*■ 0, but rather equation (11.41). 
In the limit that 3 is zero, we have 


- 63 r^n = const. (V.l) 

j 1 - 6 ,^ - 66 “ ■ 

For simplicity we assume that 3g^<<l, so that this equation 
to second order is 


1 + 




= 1 + 


1 

a-1 




(V.2) 
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From equation (11.30), 68^ in the limit 3-^0 is 


63 ^ 






(V.3) 


where Z = r/r^. From equation (III. 3) for p*/p^ , we see 
that 3^ can be expressed as a function of Z and p*/p^- Of 
course the sound velocity is 3„^ = 3„^^(p*/p*) • With 

these expressions, we can write equation (V.2) in terms of 
p* and Z alone 


eo 


ot-1 
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a-1 63. 
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(V.4) 


where is given by equation (11.51). Equation (V.4) defines 

P*/P^* a function of Z. 
o 

If a static atmosphere is to exist, we expect a top to 
that atmosphere and we now consider whether there exists a 
such that p*(Z^) = 0. If such a point exists, we see from 
equation (III . 3) and (V.4) that it must satisfy the equation 







o 


2 


(V.5) 
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This cubic equation for Z_^ does not have a real root greater 
than one if < 1. This is the regime of the thermally 
driven wind solutions. If > 0 and Xq > 1 [cf. equation 
(IV. 11 )] , there are also no real roots greater than one. This 
is the regime of the wave-driven winds , as in Section IV . 
However, if £ >0 and x '^1/ there exists a real root greater 

than one, and therefore a top to a static atmosphere. If we 
define the angle (|) such that 


90° < <p < 180° 
cos <p = - Xq^/=^ 


(V.6) 


then the top of the atmosphere occurs at 


Z 

T 


5 X 
^o'^o 


cos 


(| + 240°) 


(V.7) 


In the limit that 66^ goes to zero, Xq 90 bs to zero, and 

ij) ->■ 90° + so that is approximately («“!) ] 

This is the usual expression for the top of a gravitationally 
contained poly trope atmosphere. 

Thus, if the initial wave amplitudes are below the limit 
given by Xq = 1/ atmosphere is static, with a density 

profile given by equation (V.4) and a top Z^ given by equation 
(V.7). The static density profile plotted in Figure 2 is 
computed on the basis of these equations. In these configurations , 
the static atmosphere is modified by the presence of the waves, 
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but the waves propagate outward from r^ with no decrease in 
energy flux. They emerge from the top of the atmosphere 
traveling at the speed of light in vacuum, with 6E^ and 6H^ 
subsequently decreasing as 1/r^ [cf. equation (11.29)]. 



- 53 - 


VI. LIMITATIONS OF THE MODEL 

Under various restrictive assumptions, we have formulated 
and found solutions to a well-defined mathematical model of a 
physical situation. Before summarizing the properties of these 
solutions, we point out their probable defects and the pitfalls 
involved in applying them to real astrophysical situations. 

First, the assumption of radial fields is a strong one, partic- 
ularly in the circumstance that field pressures are large 
compared to kinetic pressures, and closed field configurations 
are to be expected. In such cases , qualitative properties of 
our model will apply only in regions where the field lines are 
close to radial for various reasons (most obviously, in the 
polar regions of a magnetic dipole) . We note that the crucial 
feature which gives rise to the possibility of wave driven winds 
is the rapid fall-off with radial distance of the density as 
compared to the field strength, and the consequent rapid increase 
of c63, the transverse perturbation velocity, with radius. This 
situation will obtain as long as H falls off inversely as a low 
power of r, and it may not be unreasonable to expect wave— driven 
expansion in the polar regions of a strong magnetic dipole. 

Short of a detailed calculation, however, we can only speculate. 

For the MHD approximation to be valid, we must have wave 
frequencies which are small compared to cyclotron frequencies 
and densities which are high enough to make the concept of a 
fluid meaningful. For the latter reason, the dynamic solutions 
just above cutoff, for which the mass flux approaches zero 
(Section IVb) , and the static solutions below cutoff, in which 
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the density goes to zero at the atmospheric top (Section V) , 
should be viewed with caution. In addition to this upper 
limit on the frequency lo, we must also satisfy a lower bound 
of the form oj >> 6V/r for all r, otherwise the transverse 
particle displacement associated with the wave will be an 
appreciable fraction of 2nr, Since 6V(r) increases much more 
rapidly than r, reaching the local escape velocity within a few 
stellar radii, this lower bound is not inconsequential, espe- 
cially for small objects with high escape velocities. 

We have ignored any possibility of wave damping, even 
though we may have a situation where the velocity per- 
turbation c66 increases dramatically with radius, approach- 
ing the escape velocity in a distance on the order of a stellar 
radius. However, we note that 6^ also increases rapidly out- 
ward, so that initially 6B/6 actually decreases with increasing 

di 

radius, as does 6H/H (see Figure 3) . In addition, there are 
reasons to believe that circularly polarized Alfven waves in a 
completely ionized, rarefied plasma are difficult to damp both 
from an observational and theoretical standpoint (Belcher and 
Davis, 1971; Barnes, 1966). 

Finally, in retrospect, one of the strongest assumptions 
we have made is that wavelengths are short compared to local 
scale heights. In many of our numerical solutions (cf. Figure 3) , 
there are regions near the top of otherwise static atmospheres 
in which wave amplitudes and radial velocities exhibit spectac- 
ular increases over short distances and the short wavelength 
approximation here becomes suspect, at the least. The situation 
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improves somewhat if we depart from an adiabatic atmosphere 
(a 5/3), as this tends to smooth out the abrupt gradients, 
and the proper inclusion of non-WKB terms would probably have 
the same effect. It is not clear how the inclusion of such terms 
would change the dynamical situation (energy and mass fluxes, 
etc.), although we would hope for no qualitative changes. To 
answer the question properly requires the numerical integration 
of the full transverse equations of motion, in conjunction with 
the radial momentum equation. 
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VII. SUMMARY 

We have investigated the properties of stellar winds in 
which the only source of energy flux is due to low frequency, 
undamped Alfven waves propagating outward along radial magnetic 
field lines. The thermal properties of the plasma are described 
by a poly trope relation, and we have considered only situations 
which would lead to static atmospheres in the absence of waves. 
We have demonstrated that Alfven waves of sufficiently large 

amplitude are capable of driving the supersonic, super-Alfv^nic 
expansion of the plasma, with a complete transfer of wave energy 
flux near the star to streaming energy flux far from the star. 

The process responsible for the acceleration of such winds 
is intrinsically different from that which produces thermally 
driven winds. To illustrate these differences, we discuss 
briefly the features of the non-relativistic solutions. Let 
us denote c6 by V, c66 by 6V, cg^ by V^, cg^ by V^, and cg^ 
t»y Vg. In the situation that V is much less than V , which 
is in turn much less than c, the differential equation for the 
radial velocity V is (cf, equation 11.48) 


r ^ 
V dr 


, V ^ - 4 (V ^ + 6v^/4) 
^ s 

^ (Vg2 + 6 vV4) - 


(VII. 1) 


In the absence of waves, this equation reduces to the familiar 
poly trope form for the radial gradient of V. In the presence 
of waves , the velocity perturbation 6V complements the local 
sound velocity V^. The dynamical effects of these two velocity 
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t 0 irnis a.ir© djco^sticslly diffciTBiit/ howsvBiT/ b 0 Cd.us 0 of thsiir 

dxff 0 x*j.n^ lD0hsviox' q.s 0 . function of jTcidisi distainc© from ths 

coron^X ]d^s0« If p xs tlis m 0 .ss donsxty ©jid ct th,© poly*** 

trope index, then is proportional to (p)^”^ . Thus V ^ is 

s 

^®st a constant (ot = 1) and at worst decreases rapidly with 
distance from the coronal base (a = 5/3) . On the other hand, 
close to the star the Alfven wave amplitudes vary so as to 
approximately conserve the wave energy flux 47rr^ V (p6V^) . 

Si 

As a consequence, 6V^ is proportional to r^ V close to the star 

(cf. oquation (11.30) for V << V^) • Sine© th© donsxty may 

decrease outward as a high inverse power of r, with H falling 

off only as 1/r^, the Alfven velocity V can increase sub- 

a 

stantially over its initial value at some reference level r • 

o 

Physically, the Alfven waves are propagating outward into an 
increasingly rarefied atmosphere, and to conserve energy flux 
initially increase outward, rather than decrease, as 

does V ^ . 
s 

The relative importance of the 6V^ and V ^ velocity terms 

s 

in equation (VII. 1) has a major influence on the nature of 
the critical point solutions of this equation. In the absence 
of waves, Parker (1963) has shown that equation (VII. 1) will 

have well-behaved wind solutions for initial values at the refer- 
ence level r^ in the range 


a - 1 


V ^ 
so 


eo 


< 4 


(VII. 2) 


The lower limit in this inequality represents the point below. 
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which critical point solutions to equation (VII. 1) no longer 

exist (the gravitational field is too strong to allow expansion 

and the corona assumes a static configuration) . The upper limit 

represents the point above which conditions at the coronal base 

are no longer realistic (e.g. , the streaming velocity is already 

supersonic at r^) . Consider the situation in which the lower 

limit in inequality (VII. 2) is violated. In the absence of waves 

there will be no critical point solutions to equation (VII. 1) 

because V^^ is never comparable to . With waves, however, 

there is still a possibility of a wind solution if 6V^ at the 

critical point is comparable to V^^ there. Let the subscript 

"c" indicate evaluation at the critical point r (if it exists) . 

c 

The above, discussion of the radial dependence of 6V^ implies 

that for 6V ^ to be comparable to V ^ , we must have 
C 0 c 

^^o^ ^^c^^o^ ^ ^^ac^^ao^ comparable to V^^^ ’ Critical 

point solutions of this nature do in fact exist, and usually 
occur within a few stellar radii. For a given 5V^, we require 
that the density at r^ be low enough (and thus the value of 
V /V high enough) so that 6V is comparable to V . As 6V 

oG aG C 0C O 

becomes smaller, the density profiles become more rarefied, so 

that 6V still attains the escape velocity at r . However, 

c 

this process cannot continue indefinitely, since the maximum 

value possible for 6V ^ is on the order 6V ^ (c/V ) , corre- 

c o ao 

spending to the limit that V is close to the speed of light 

cl C 

(and assuming that r is on the order of r ) . If SV ^ (c/V ) 

o o o 

is substantially below V^^^ , 6V^ can not reach the escape 
velocity at r^, and critical point solutions do not exist. 
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We note that Alazraki and Couturier (1971) and Belcher (1971) 
allow arbitrarily large Alfven velocities when considering, 
the problem. In such a situation, 6V can reach V no matter 
how small <SV^, as long as it is non-zero. As a result, these 
authors incorrectly conclude that wind solutions are always 
possible if > 0. 

These heuristic arguments provide some qualitative insight 

into the range of initial wave amplitudes necessary for the 

production of Alfven winds. If we assume that V and V 

eo ao 

are small compared to c, then our detailed mathematical 
treatment indicates that well-behaved wind solutions exist for 
initial wave amplitudes in the range 

.V 6V ^ V 

It ”1^ < < MIN(1, (VII. 3) 

^ ° ^ V ^ ''eo 

eo 

where 4 is a factor of order xinity for our purposes [cf. equa- 
o 

tion (11.51)]. We emphasize that inequality (VII. 3) only applies 
when the lower limit in inequality (VII. 2) is violated (e.g., when 
there are no Parker wind solutions). Under various assumptions, 
we have obtained approximate analytic solutions to the critical 
point equations for initial wave amplitudes in the range defined by 

inequality (VII. 3) . We have referred to this range (which 
does not include the lower limit) as the intermediate airiplitude 
range. Initial wave amplitudes just above and inclusive of the 
lower limit of inequality (VII. 3) are said to be in the cut- 
off regime, as discussed below. Initial wave amplitudes below 
this lower limit are too weak to drive the coronal expansion. 



-60- 


and the atmosphere assumes a gravitationally bound, static 
configuration. Wind solutions exist for initial wave ampli- 
tudes above the upper limit in inequality (VII. 3) (referred to 
as the strong amplitude regime) , but they are such that either 

V > V (if V > V ) or V > V (if V < ) . The 

o ~ eo ao eo c - ac ao eo 

situation in which V > V is clearly unrealistic for a 

o ~ eo 

reference level at a coronal base. Although the situation in 
which V > V at r does not necessarily imply that conditions 
at r^ are unrealistic, it does prevent us from obtaining 
approximate analytic solutions to the critical point solutions. 
We must rely on numerical solutions for wave amplitudes above 
this limit, and as a consequence our understanding of the 
properties of solutions in the strong amplitude regime is 
limited. For wave amplitudes in the intermediate and cutoff 
regimes, however, our understanding is detailed - 

Before considering these properties, we note as cin aside 
that our choice of parameters for the specification of the 
initial conditions at r^ is neither unique nor necessarily 
ideal. For example, since 6V and V vary rapidly as functions 
of radius , it may be difficult in practice to choose values 
for them at a given point, since they are sensitive functions 


of distance from the star. A quantity which may be more 


appropriate as an initial value is the Alfven wave Poynting 

flux at r , 4Trr ^ (c 6E 6H_^/4 it), since this flux is essentially 

constant close to the star. If we denote this flux by t 

let W be 4irr ^ V (H ^/8 tt) , then over the range of wave 
o o eo o ' ^ 
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amplitudes in equation (VII. 3), F ° varies between the limits 

P 


^ £ 
2 7 


V 

3 eo 


W 


V. 

< 2W^ MIN(1, ^ 


eo, 


(VII. 4) 


ao 


Well-behaved wind solutions will exist if the Poynting flux 
associated with the Alfven waves at lies in the range given 
by expression (VII. 4). This requirement on the Poynting flux 
for the existence of well-behaved wind solutions is completely 
equivalent to the requirement on the wave amplitudes given 
above . 

The properties of stellar winds which are of primary 
interest are their energy fluxes, mass fluxes, and energies 
per particle at infinity. For initial wave amplitudes in the 
intermediate range, the streaming energy flux at infinity is 
equal to the Poynting energy flux at the reference level, 

In this range, the energy flux at infinity therefore varies 
between the limits of expression (VII. 4) and is proportional 
to 6V^ . The mass flux F^^ in the intermediate range varies 
between the limits yiven in expression (IV. 34) and is proportional 


6V 


The lower limit in this expression does not include the 


abrupt decrease in mass flux as 5V^^ approaches cutoff. The 
limiting velocity y„ of the wind for the intermediate range 
varies between the limits given in inequality (IV. 36) and is 
proportional to 1/6V^ . Again, the upper limit in this equa- 
tion does not include the violent behavior in the regime of 
wave amplitudes just above cutoff. The distance required for 
a significant transfer of the initial wave energy flux F 


o 
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to wind streatming energy flux is of the order of r^, where 
r^ is the point at which the radial streaming velocity is 
equal to the local Alfven velocity (Section IVe) . For the 
intermediate range of wave amplitudes, assuming 
r varies between the limits given in inequality (IVo45), and 

di 

— V 

is proportional to 6V^ In the intermediate range, 75% 

of the kinetic energy per particle at r is associated with 
the radial streaming velocity, and 25% is associated with the 
transverse velocity of the wave. 

As closely approaches the lower limit in expression 

(VII. 3), the mass flux decreases abruptly toward zero, and 
since the energy flux remains well-behaved, the energy per 
particle at infinity increases abruptly toward infinity. In 
the same limit, the distance r goes to infinity. Qualitatively, 

a 

this situation arises from the fact that just above cutoff the 
Alfven waves are able to drive a vanishingly small mass flux off 
the star. At cutoff, we have a finite wave energy flux at r^ 
which is to be distributed at infinity among an infinitesimal 
number of particles. Below cutoff, dynamic expansion is no 
longer possible, and the atmosphere formally assumes a static 
configuration. In the static case, the energy flux associated 
with the waves is rigorously conserved, and appears at infinity 
in the form of waves. As pointed out above, the physical sig- 
nificance of the dynamic solutions just above cutoff, and the 
static solutions below cutoff, are questionable because of the 
extremely low densities involved. 
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For the purposes of nximerical illustration, consider 
a star with V ^/V ^ = 0.1, a = 5/3, r = 0.75 x lo ^ ^ cm, 

= 599 km/sec, and = 0.5 gauss. From expression (VII. 4), 
the initial Alfven wave Poynting flux needed for the pro- 
duction of an Alfven wind lies between the limits 0.85 x 
10^** ergs/sec and 0.84 x ergs/sec. If we take a density 

at r^ of 1 X 10* particles/cc , the Alfven velocity is 

109 km/sec, and the intermediate range of wave amplitudes varies 
between 2.3 km/sec and 256 km/sec. For wave amplitudes in 
this range, we obtain the same streaming energy fluxes at in- 
finity as Poynting fluxes ar r^, and mass fluxes between 
the limits 3.1 x lO® gms/sec and 2.9 x 10^** gms/sec. The 
streaming velocities at infinity lie between 24,500 km/sec 
and*^I3ip km/sec. 

In conclusion, we note several points in comparing Alfven 
winds with thermally driven winds. To obtain winds in the 
thermal polytrope models, the sound velocity at the coronal base 
must be on the order of the escape velocity there, and the 
limiting velocity at infinity is of the same order. In Alfven 
winds, the initial transverse velocity 6V^ can be small compared 
to the escape velocity, by as much as a factor of order v'V^ ”/ c , 
and the limiting velocity at infinity can be large compared to 
the escape velocity, by as much as a factor of order /o/V^^ . 

Even in tightly bound atmospheres with small initial wave ampli- 
tudes (6V ^ << V ^ << V ^), rarefied, energetic wind solutions 
o so eo 

may exist. The existence of such winds even in these extreme 
situations is closely related to the collective nature of the 
acceleration process, in which the energy input required to 
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maintain the low velocity, transverse motions of a great many 
ionized particles at a coronal base is ultimately transferred 
to the high velocity, radial streaming of relatively few 
particles far from the star. As a result, wave driven winds 
can exist even in tightly bound situations, and may exhibit 
relatively high energies per particle at infinity. It is this 
qualitative concept of the acceleration mechanism that we wish 
to emphasize, rather than the quantitative details of the solu- 
tions we have presented. 
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Appendix A 

We solve equations (11.22) and (11.24) for the WKB wave 
amplitudes under the assumptions: 1) the wavelengths are 

small compared to the local scale heights; 2) 63 and 3g are 
small compared to one, so that we may neglect second order and 
higher terms in 63 and 3^. With these assumptions, we can 
immediately drop the last term on the left hand side of 
equation (11.24), and replace F everywhere by y# where 


Y = 


(n - ^) ’^ 


(A.l) 


Strictly speaking, to justify this approximation, we must re- 
quire that 63 be small compared to 1/y rather than 1, but for 
our purposes y can be taken to be close to unity. 


We may also replace (e + p) in equation (11.24) by p*c^ , so 
that equation (11.24) becomes 

[y^ I It If 

(A. 2) 


Hr , ^ 
Tir ^ cri 


SE ifj (r6H)l 


0 


We eliminate 63 from this equation by using equation (11.21) 
in the form 

40 = - SE 1 (A.3) 

r 

For convenicnco, we introduce a new variable 6h defined 

by 

6h = ^ (A.4) 

n 
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Using equations (A. 3) and (A. 4) in (A. 2), we obtain after some 
raanipulation an equation for 6 h and 6 E alone, as follows 


H 


[ Y' 


4lTp*C^11 


] 11^ 6E + 2y’6 ^ 6E 


+ 


nH^ 

r 

47T p*c^ 




] |- 6 h = - 6 E( 6 y 



4By^ 

r 


(A. 5) 


+ 6h [ 


dr 


H 


47rp*c^r 


3y^ B^ 
r 


+ y3^ 



] 


Equation (11.22) is our second equation for 6 h and 6 E, and 
has the form 


1 3 ih + I- «E = - 

c 9t dr r 


(A. 6 ) 


We now apply the method described by Weinberg (1962) 
to obtain the WKB amplitudes (see also Belcher 1971b). Let L 
be the scale height for variations in p*, B, and H^. We write 
5E and 6 h in the form 


6 E(r,t) 
6h (r , t) 

where 


[ 6 Ej^(r) + p 6 E 2 (r) + p^ 6 E 2 (r) + ...)exp[i(wt - S (r) ] 

[ 6 h, (r) + y 6 h-(r) + y^ 6 h-(r) + ...]exp[i(wt - S(r)] 

(A. 7) 
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dS 

dr 
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2v 

kL 


(A. 8 ) 
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The parameter y is the ratio of the wavelength to the scale 
height L, and is assumed to be small. The quantities '^^2' 

iShj^, etc., are also assumed to have scale heights on the order 
of L. We insert equation (A. 7) into equations (A. 5) and (A. 6) , 
and keep only terms to first order in y . If we place zeroth 
order terms in y on the left hand side, and first order terms 
on the right hand side, we obtain the equations 


_ |li» 6hj^ + ik = - y ^ Sh^ + ik SE^] 


lOJ 


6E 


1 d 


j <SE, 

r dr 1 


(A. 9) 


and 


H 


- iiiJl (Y^ + ) 6E, + i2y^Bk 6E. 

4Trp*c^n 


nH ^ . H ^ 

+ ik ( £ 32 y"") <Sh, = - y[- ^ (Y^ + ~~) 6E2 

4iTp*c^ 4irp*c^Ti 


nH, 


+ i2y^ek 6E„ + ik ( S^Y^) 5H-] 

4ttp*c* 


+ 4h^ tSY“ ^ - 
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4irp*c^r 




- 6E^ [ST a? 


HH 


- 3F «*■! 


4ttp*c 


{A. 10) 
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where we have assumed that r > L. The zeroth-order approxi- 
mation is obtained by assuming the right hand sides of (A. 9) 
and (A. 10) are zero. Using 6 as defined by equation (11.26), 

3 . 

we find that 


3 


P 



0 + 110 ^ 

1 + 33a/n 


(A. 11) 


and 


6 h 


1 



(A. 12) 


We have chosen the sign in equation (A. 11) for outwardly 
propagating waves. 

To obtain the first-order solutions, we insert the zeroth 

order solutions (A. 11) and (A. 12) into equations (A. 9) and 

(A, 10). To eliminate the quantities 6 E 2 and <Sh 2 , we then 

multiply equation (A. 9) by (nH V4irp*c^ - and add it 

P 

to equation (A. 10). This leaves us with a differential equation 
for (SEj^ of the form 


57 ‘ 7777 - i? ^ a? 

3 


+ { 


1-3. 


- r 2 . 




6 E 


i (3V^ g 


1 - 6 , 


_ i + + ^2 ax, , 0 

2 r r ' dr 


(A. 13) 





-A5- 


To solve equation (A. 13), we write all derivatives cf the 
form dp/dr, dy/dr, dg^/dr, and so on, in terms of dp*/dr. For 
example, we have from equation (A.l) that 

r & (A.14) 

dr ^ I n dr 2 ' dr ^ 

We also have [cf. equation (11.18)] 


2 a. 1 + i ^ ^ = 

r p* dr Y dr 8 dr 


0 


(A. 15) 


From (A. 14) and (A.15) we can obtain expressions for d6/dr 

and dy/dr in terms of dp*/dr alone* From the definition of 

3 , we also can obtain an expression for d^^/dr in terms of 
a 

dp*/dr [cf. equation (11.43)1, Proceeding in this manner, we 
eliminate all derivatives in equation (A. 13) except those of 
p* and 5Ej^. After a tedious process, we obtain the equation 




dp* 

dr 



(A. 16) 


The solution to this equation is (11.29). We obtain expressions 
for the radial dependence of and 68^^ by using equation 

(11.29) in conjunction with equations (A. 3) , (A. 4) and (A. 12). 
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Appendix B 


We sketch a derivation of the non-relativistic equation 
(11.48) iri the limit that $ and are small compared to one, 
with no restrictions on 3^- The possibility that 3^^ may be 
close to one means that we must keep time derivatives of E 
in Maxwell’s equations. For a non— relativistic MHD plasma 
with fluid velocity V, density p, magnetic field H, and 
electric field E, in the presence of a spherically symmetric 
gravitational potential 4*, the relevant equations are 


pgt Y + :?p + 


pV4> 


+ ^ H 

4it ~ 


[V X H - 




= 0 


(B.l) 


V X E 


+ - -Ir H 

C dt - 


= 0 


(B.2) 


E + — V X H = 0 (B. 3) 

-- c 

If we assume that V = V(r)e^ + 6V{r,t)§^ and 

H = H(r)§^ + 6H(r,t)§^ , then from equation (B.3) we immediately 
have E = 6E(r,t)e_, with 6E = (V6H - H6V)/c. The 4>~component 
of equation (B.l) is 


9 

9t 


6V + 


V 

r 9r 


(r6V) 


_H ^ 

4iTpr 9r 


_H ^ 

47TpC 9t 


(r6H) + 


6E 


(B.4) 
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The 4>-component of equation (B.2) is 

Thus , our perturbation quantities satisfy equations which are 
the non-relativistic limits of equations (11.21), (11.22), 

and (11.24). Hence, the WKB solutions to equations (B.4) and 
(B.5) are the non-relativistic limits of the solutions given 
by equations (11.27) through (11.32). For convenience, we 
now assume our perturbations are circularly polarized in the 
equatorial plane of a spherical polar coordinate system. 

The radial component of equation (B.l) may be written 
as 


V ^ - ixi 4. 1 ^ ^ 

dr ~ r p dr dr 4ttpc 3t 


Ir (E 


X H) 


(B.6) 


1 

4iTp 


{H X (V X H) + E X (7 X E) = 


where we have used equation (B.2). Our assumption of circular 
polarization implies that E x h is time-independent, so that 
equation (B.6) becomes 


V 


dV 

dr 


1 ^ 
p dr 


dr 


rp 


l~ (iSH^ + 6E^) - pSV^] 


(B.7) 


+ 


1 

8ttp dr 


t6H^ + 6E^3 


= 0 


- 7 /- 
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We may write equation (B.7) in terms of dS/dr alone by using 
the conservation of mass and equation (11.29) for the radial 
dependences of the transverse perturbations in terms of p, V, 
and r. In particular, we note that in this limit 



(1 - 3 ^") 

C2. 

(1 + 83 ^^" 


dr 


(28 + 3^) - 



(B.8) 


After some tedious manipulations, we recover equation (11.48) 
for dV/dr, neglecting third order terms in small quantities. 
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FIGURE CAPTIONS 



Figure 1. Critical point solutions (r /r , g /3 / 3 /B 

^ ^ ' c o c eo o eo 

as functions of values of ^gQ^/^eo^ 

equal to 0.1, 0.040, and 0.001 (curves labeled A, B, 

and C, respectively). The Alfven velocity 

and the escape velocity are constant, with 

3^^ = 0.0001 and = 0.0020. The limiting velocity 

of the wind at infinity is cB^ . The dotted curves 
represent analytic results derived in Section IV 
under various approximations (see text) . The scale 
for r /r^ is to the right and the scale for all 
other quantities is to the left. 


Figure 2. Energy flux and mass flux as functions of 

for the same values of ^gQ^/^eo^ ' ^ao' ^eo 

in Figure 1. We also plot F^/c^Fj^ - 1, which is 

y - 1. 

Figure 3. The curves labeled 1 through 4 correspond to the 

values of 6B ^/B ^ indicated on curve A of B_/B_ 
in Figure 1. The curves labeled 5 are for a value 
of indicated by a vertical stub on the axis 

in Figure 1 (this point is below the cutoff for dy- 
namic solutions). The quantity B ^/B^ is 0.1, with 

so eo ' 

6 = 0.0001 and B = 0.0020. We plot different 

ao eo 

variables as functions of radial distance from the 





reference level r^: (a) The mass density normalized 

to its value at r^, and the radial streaming velocity 
divided by the escape velocity at r^; (b) The Alfven 

velocity divided by c, and the transverse velocity 
divided by the Alfven velocity; (c) The transverse 
velocity divided by the escape velocity at r^, and 
the transverse velocity divided by the radial velocity, 
(d) The Alfven wave energy flux A normalized to its 
value at r^, and the radial Alfvenic mach number. 

The scales for the various quantities are to the 
left or to the right, as indicated. Vertical stubs 
on the curves indicate the location of the critical 
points . 
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FIGURE 3a 












